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ABSTRACT 


This paper investigates the application of the Kalman Filter and the 
General Exponential Smoothing techniques of forecasting. Both methods 
are derived and the similarities and differences between them are dis- 
cussed. The two techniques are then applied to the practical problem of 
predicting weekly losses suffered by the U. S. Marine Corps units in the 
I Corps Tactical Zone in the Republic of Vietnam. The mean absolute 
error of the prediction is used as the criterion for choosing the better 
of the two methods. Results are given for both techniques as well as 
for the method of linear regression. In general the Kalman Filter pro- 
vides the smallest mean absolute error for the three mathematical models; 
linear, growing sine with harmonics and frequency of sixteen, thirty-two, 


and fifty-two weeks, and a constant model. 


Library 


U.S. Naval Postgraduate SchooT 
Monterey, California 93940 


TABLE OF CONTENTS 


= INTRODUCTION --------------------------------------------- 7 
A. THE MODELS --------=-=============----====----------- 8 
B. NOTATION --------------------------------------------- 8 
II. THE KALMAN FILTER ---------------------------------------- 9 
A. THE KALMAN FILTER EQUATIONS -------------------------- 1] 
B. REMARKS ---------------------------------------------- 12 
III. GENERAL EXPONENTIAL SMOOTHING ---------------------------- 14 
A. THE WEIGHTING VECTOR --------------------------------- 15 
B. THE FORECAST BIAS ------------------------------------ n 


C. EXPONENTIAL SMOOTHING COMPARED WITH KALMAN FILTER ---- 22 


ieee SPORECASI IN GENER INGMROSSES 5222222 e 23 
A. FORECASTING WITH EXPONENTIAL SMOOTHING --------------- 24 
BT rORECAST ING WITH THE KALMAN PIETER ER 42 
С. FORECASTING WITHMIEINEAR REGRESSION ----------- ee 44 
EROGO See Em Жез nn ee 69 
A A ee 74 
O O: o σε EE 79 
O A ο Ш 80 
O E a ee ee eee eee — 8] 








ων, 


III. 


Vale 


VT 


ЕШ 


De 


Ale 


LIST OF TABLES 
GENERAL EXPONENTIAL SMOOTHING RESULTS, GROWING SINE 
WITH HARMONICS, CYCLE LENGTH 16 WEEKS ------------------ 


GENERAL EXPONENTIAL SMOOTHING RESULTS, GROWING SINE 
WITH HARMONICS, CYCLE LENGTH 32 WEEKS ------------------ 


GENERAL EXPONENTIAL SMOOTHING RESULTS, GROWING SINE 
WITH HARMONICS, CYCLE LENGTH 52 WEEKS ------------------ 


GENERAL EXPONENTIAL SMOOTHING RESULTS, LINEAR MODEL ---- 
GENERAL EXPONENTIAL SMOOTHING RESULTS, CONSTANT MODEL -- 


KALMAN FILTER RESULTS, GROWING SINE WITH HARMONICS, 
CYCLE LENGTH 16 WEEKS ---------------------------------- 


KALMAN FILTER RESULTS, GROWING SINE WITH HARMONICS, 
CYCLE LENGTH 32 WEEKS ---------------------------------- 


KALMAN FILTER RESULTS, GROWING SINE WITH HARMONICS, 
CYCLE LENGTH 52 WEEKS ---------------------------------- 


KALMAN FILTER RESULTS, LINEAR MODEL -------------------- 
KALMAN FILTER RESULTS, CONSTANT MODEL ------------------ 
LINERRRECHESSTONMRESUETS = ERE NN 


сл 


27 


30 


33 
36 
29 


45 


48 


21 
54 
57 
60 





I. INTRODUCTION 


One of the problems in any forecasting technique is what to do about 
additional data. That is, given a forecast based upon some data, how 
does one go about updating this forecast when additional data is obtained? 
Two such methods are discussed in some detail. These are the methods 
known as the Kalman Filter, named for R. E. Kalman, and General Exponen- 
tial Smoothing, a name coined by Robert G. Brown. 

Both the Kalman Filter and General Exponential Smoothing are derived 
and comparisons made for any similarities and differences. Both tech- 
niques are then applied to the problem of forecasting losses incurred by 
the U. S. Marine Corps operating in the I Corps Tactical Zone of the 
Republic of Vietnam. The actual losses are expressed in weeks and cover 
the period from December 24, 1965, to April 6, 1968. Losses include al] 
Marine Corps personnel who die as a result of battle as well as non-battle 
causes. It also includes those personnel who must be evacuated from 
Vietnam because of wounds or injuries received from any cause. 

The basis for comparing the Kalman Filter and General Exponential 
Smoothing lies in their similarities. In each case, a mathematical model 
1s assumed. The model consists of coefficients, which are unkown, and 
functions, called fitting functions, whose values are known. The problem 
becomes one of estimating the coefficients, based upon previous data, 
then using this estimate with the fitting functions to obtain a forecast. 
The forecast error, the difference between what actually occurs and what 
has been forecasted to occur, is weighed to obtain a new estimate of the 


coefficients. 


A. THE MODELS 
The models used for forecasting U. S. Marine losses are a linear 


model, 


PERE a, + apt, 


a growing sinusoidal with harmonics, 


2 Ont 


С - (а, + а Е) + (a, + a,t)sin ^ 





¿mt "mt Ant 
(ад + а Е)со5--- + a, sin 4 42 СӨ5----, 


and a constant model, 


In the second model, the forecasts were made for values of w of sixteen, 
thirty-two, and fifty-two weeks. 

The second model was chosen since there appeared to be some cycle in 
the number of losses incurred. It was felt that losses would increase 
as a result of increased enemy activity. Such activity was at a low ebb 
during the time of the monsoon rains. During the dry period, the tempo 
of enemy offensive actions increased for a time with a resultant increase 


in Marine Corps losses. 


B. NOTATION 


The notation used above reflects that used throughout this paper. 


The symbols a(t) or C(t) represent the forecast of losses for the pth 


time period. The same symbols without the above them are the observed 


value or the actual losses during the tth 


period. It is also convenient 
to adopt vector and matrix notation. Thus, an alphabetical symbol, such 


| Ж 
as a, will represent a column vector, while a will denote a row vector. 


In general, a vector with a subscript will represent the vector at a 
Specific point in time. Matrices will be denoted by a capital alpha- 
betical character. In the event the matrix symbol is subscripted, it 

will denote the matrix obtained by taking the expected value of a col- 
umn vector, the first subscript, multiplied by the transpose of the vector 
denoted by the second subscript. A matrix with a single subscript is the 
expected value of the vector represented by the subscript times its trans- 


pose. For example, 


ab 


II. THE KALMAN FILTER 


Basic to every forecasting technique is the criterion which is used 
to determine the forecast. In the case of the Kalman Filter, which esti- 
mates the coefficients of the model, the criterion is to minimize the 
trace of the estimation error covariance matrix, Ca where e is the 
unobservable error vector between the true value of the model coeffi- 
cients and the current estimate of these coefficients. This estimate is 
the minimum mean square error estimator. 

A fundamental theorem of estimation, called the Gauss-Markov Theorem, 
is very important in any form of estimation. The theorem shall be stated 
but no proof is provided. 

Theorem: If o is the vector of the observed data at times t. and x 
is a linear function of o, X and 9 are random variables with moment matrices 


- ] | Р E 
C. and C , and if C exists, then the linear minimum mean square 


7 Χρ 


estimate x of x given ) 15 


Cast G 


and (е X ХО 779 


Throughout this paper, the data will be assumed to be linearly 
related to the vector x (or a). Thus, 


gå =) Dx +, 


— 


where v is the measurement error or noise vector and B is anm x n 
matrix whose rows are the vectors είς), the vector of fitting functions 
evaluated at time t. For example, if o(t) = X} + X,t + v, the vector 


of fitting functions f(t) = (1,t) and the vector of coefficients 


l: 
w (x, я хо). 
Then, 
9 = E[x(Bx + ЯШ 
хө ur 
» T 

as C o & СР t бху 

_ Т 
апа Саз == В te Breie VINN: 

ai T ПЕ Т 

or Co? = ВС, В ar BC, + CB + С, я 


From the fundamental theorem, 


Т Т -1 


. _ T 
¿sssi [CB + C ПВС,8 + Весе + б, + CJ O > 
I Τ EL т, T 
and C sinis С, - [C 8 + С ЈС8С, В + TE + BC, + CJ [С, В + б] Е 


Normally the vectors x and v are assumed to be independent. This, 
with the assumption that E(v) = 0, implies that C uis is identically zero. 


With this assumption, the above equations reduce to 
| а) ТА . 
X C.B [BC_B + Cy) 


10 


λα. 


_ T T -]γρρ]η] 
and C = C, - СВ [BC_B di С] [BC] 


e 


By the use of the matrix identities 


21 1.-i 


(2.1) АВПС > BABI]! = A”! + Bicla” gle 


TT (2.2) [C «la !gr! s c cal ов! 1 ac, 


the following equations are obtained: 


а ie A B'c-'g] "ВГС Т ole 


g -] Ta-1p3-1 
and Co = κ + В С, Bde . 


There still remains the matrix C, which must be known before the 
estimate x can be made. In general little information is available 
concerning x. In this case it is reasonable to assume that the diagonal 
elements of ο are larger than any finite number. With this assumption, 


an is considered to be 0. Then, 


- [B ву 1 ВГС o 


= 
II 


Τ--]ρη-] 
[Β С, В] 


and 


A. THE KALMAN FILTER EQUATIONS 

The foregoing has been the groundwork. Suppose that m observations, 
o", have been taken and т is the minimum mean square error estimate 
based upon the m observations. The matrix = is available. In addition, 
suppose that r additional observations are taken. The problem is to 
determine the estimate ІЛЕ 
Since и is available and is 


^ 


M = ΠΠ gm! ВИТ (С ο. Kram 


eu E rg? T ( cm 16117-1 


e V 
m m m 
Also, 0 = Вх + ү е 
апа mer mtr mtr πι. 
Β X T 


^ ^ ^ 


The noise vector v^ is independent of the vector v. This implies that 


СШ is a diagonal matrix with elements e and ος, By partitioning the 


er 


^Y. 
matrix BY! into components B" and B” and the vector o" into 6" and o 


the new estimate 


à Τ -1 
omr = К" E BT A: г” κ ixm + ΒΒ (ου) 9" 


is obtained. Also, 


mtr 4 (сту! Å g^ T (С) lgrz7! | 


6 е V 
Then define K = qum . From equation (2.1), 
e cr + pre" И. 


From equation (2.2), 


Finally, the new estimate is 


- қ : 
es Ща" k[o* ` B 


^ ^» ^ = ^ 


These last three equations are the Kalman Filter. [Liebelt 1967, 


"mtr 


pp. 165-166] Note that succeeding estimates x are computed recursively 


қ > ^m 
from the prior estimate x . 


B. REMARKS 

Notice that with the assumption that С, is known, neither K nor Ce 
depend unon the data 0. These elements could be precomputed and re- 
tained for later use as needed. Also, it is not necessary to assume 
that the matrix c. is zero. If one has information about the distri- 
bution of X, the elements of C. can be computed and used in determining 
"E 

The vector K may be thouaht of as a weighting factor applied to the 

forecast error. Notice that the weight for any error is proportional 
to the inverse of the variance of the forecast. A large variance in 
the measurement or noise for any particular forecast will cause the 
error to be given less weight in determining the estimate of the coef- 
ficients to be used in computing the next forecast. 

The dimension of the matrix б) is always r x r. If measurements 
or observations are taken at each succeeding time period, сү is a scalar 
quantity and is the variance of the measurement for that period. In the 
event that the matrix С, is diagonal with the diagonal elements equal to 
some constant raised to some power, the matrix Ce will reach a steady 
state. That is, 


mer m 
= C 
Ce е 


The power to which the constant is raised corresponds to the time period 


for which the observation is made. 


III. GENERAL EXPONENTIAL SMOOTHING 


As is the case in the Kalman Filter method, the linear model 


is assumed, but rather than estimating x in order to minimize the 
diagonal elements of some matrix, it is desired to minimize the sum of 


the weighted squared error, 
m 
2 i КЕ 2 
> w.[6(j) - 6(j)] . 
j=l 
If W is a diagonal matrix with elements (Wise ит), then 


m - på = 
şi νἴ[ο(1) - (1) ]2 = [We - 6) Ce 1. 
i= μα. ma 


Also, it is assumed that 


where the matrix B is defined the same as in the discussion of the Kalman 


Filter. 


Define the matrix F to be ВВ. Then 


ts Fe glass is al 


I 
- wer ro ва + ШИ 


The matrix F is positive definite. Therefore the choice of д! which 


minimizes the above expression 15 


АТ = wer I 


Finally, 


ама με 


= 


1 


Notice that this is identical with the Gauss-Markov estimate with 


_ qe Т ¿3 
ЕВА, =o 


B and WW 
A. THE WEIGHTING VECTOR 

A fundamental hypothesis of Exponential Smoothing is that the 
vector of fitting functions {({11) must be a linear combination of f(t). 
The matrix of coefficients of this combination, which does not depend 


on time, is called the transition matrix, L. 


From this, 
πο = EFE 

and tem) = FI(1)LIx(t-1) , 
A Е Т t 

where KE) = O 


-- ^ 


The fittina functions are always computed at the number of time periods 
in the future for which the forecast is desired. This implies that the 
Dresent period is always considered as the period t = 0. Then the 


problem becomes how to compute x(t) recursively from x(t-1) so that 


b'[o(t-i) - x(t) F(-1) 12 is minimized. 


и ct 


120 


Define the matrix 


T 
F(T) = 7 
i=0 


and the matrix F to be the limit of F(t) as t becomes very larae. The 


b'r(-i)f (3) , 


matrix F will exist if each fitting function f. (t) is areater than the 


square root of Do 
15 


Let the data vector be 


| er 
g(T) = у БОТ * 
т і-0 * 
th | TE 
Then the i^ element of the vector g(T) is Ь e(T-k)f ; (-k) 
k=0 
Then g(T) = 9(T)f(0) + bL7'g(T-1) , 
^ T 1 di ° 
and ipl) i af akai T) 
Е 1Ξ0 Е τ Я 
ЕСИ. 
But o(T-i) = f (-i)x(T) 
Hence, F(T)x(T) = oer 


Substituting the steady state matrix F and combining the two equation 


for g(T), the expression 


X(T) = (ТЕ f(0) + bF' | L lFx(T-1) 
is obtained. [Brown 1963]. 
But pF LTTE = LU -F rolf (ОЕ 
and бі) = flo. 


1 


Define the vector h to be F f(0). Then the expression for x(T) becomes 


^ 


mm > FTV. 


or X(T) = L'x(T-1) + h[e(T) - (T). 


This last equation is the Exponential Smoothing means of estimating the 


model coefficients. [Brown 1963, pp. 174-177]. The vector h is the weight 


of the forecast error applied in obtaining the new estimate x(T). Note 


that x(T) is computed in a recursive way from x(T-1). 


16 


Β. THE FORECAST BIAS 

One may question the accuracy of the forecast obtained by means of 
Exponential Smoothing. It has been shown [Bessler and Zehna 1968] that 
for the model Cit) = ar ait) , 
where E[C(t)] = a 


that the forecast is biased în the early stages by a factor of E 

where b is the weighting scalar applied to past data. Since b is 

between zero and one, then E[C(t)] = a 

when bt has vanished. However, this shows that even where the mathe- 

matical model is correct, there is an initial bias in the forecast. 
This may lead one to ponder what occurs when the assumed model is 

incorrect. It has been demonstrated [Brown 1963, p. 128] that in the 


case where the assumed model is 
E[C(t)] = a 
but the data is actually corresponding to the model 


Ecc = a, tat, 


that E[C(t)] = аре ant = a,b(1-b)"' å 


after pt vanishes. Thus, even after a sufficient amount of time has 


elapsed so that pt = 0, the forecast is incorrect or biased by a constant 


times the slope of the line through the data. 
The next step was to determine if these results could be generalized. 


To do this, the Exponential Smoothing equation 
mE Tena 
a(t) = L a(t-1) + h[C(t) - f (1)a(t-1)] 


was examined. This can be rewritten as 


17 


a(t) = iz [Uo hf he) * D - nr 0730). 
Then if E[c(t)] = gi(t)a, 
х t-] 5 
E[C(t)] = gene + АТ ЕТ-ВР (1) [а(0)}. 


The above expression is true in general when the assumed model is 


c(t) = fi(t)a + e(t) 


~ 


and the data actually follows the model 
СЕ quia ae t). 


In order to compare this expression with those obtained previously for 


the constant model, the values 


——— 
tle 
and g(t) = 1 


are substituted to obtain 


E t-l 
μον чу 
k=0 


Пн на. 


This can be written as 
Б Е L 
E[C(t)] » a[1-b'] , 


where b is I-h. Thus the result is the same as before. 
Next, it was desired to consider the situation where the data was 


modeled by E[c(t)] = a , but the data actually obeyed the model 


18 


E[C(t) = a, + aot. 


Then substituting the values 


EA 
fU + Г. 
91%) = (1,8), 
and I a! = (а, а.) 


into the expression for E[C(t)], 


К t-| k r 
ο st) ат в. 


This becomes E[C(t)] = (a *ast) (1-h*) - à СЕГИС | 


1 д 


When nt vanishes, the expression becomes 


1 


E[C(t)] = а + apt - ap(I-h)h”, 


2 2 


which is identical with the prior result with b = ]-h . 

Notice that the expression for E[C(t)] has been used to determine if 
a bias exists when the correct mathematical model is assumed as well as 
to determine the amount of the bias when an incorrect model is assumed. 
It was desirable to determine if the forecast was biased when the assumed 


linear model was correct. The equation then becomes 


ЕССЕ) = (aytayt)# (1)(I-L'+hf! (1) [1-(LT-hf!(1)) "Th 


2 


- afe) pmec er (1) ағ (АП -А ағ (191 


(1- (LI -ne 0) в Ge) Gc nre qo rr! i-r Е | 


It became important to know if the vector me. vanished as 
t became large. If the vector fails to vanish, the forecast will remain 
biased. The vector will vanish if the eigenvalues of Lene! are 
all between minus one and plus one. То demonstrate this, let 


А = Denn where A is n x n and nonsingular. Also, let the eigen- 


values of A be ШЕ Mos + > „Мо where each m. lies between minus one 
and plus one. Let Мү Ма» н У be eigenvectors of A. These vectors 
form a basis for pu 


Then AV, = M Vy s 


-- 


Also, MN Ela then 
X = ату т ...® а у 
Premultiplying each side by A, 
Ах = а] АУ] TCU + а АМД: 
Substitute to obtain 


x = amea F í m v. 
Ах ] 1.1 n n.n 


Repeating this process t times, 


t 


Since this is true for all x in a Αη = 0. 


20 


In the case of the linear model, the eigenvalues are equal to 


-] 
1 - (h. + h,)2 . 


one if h- + h; lies between zero and four. The h for the linear model 


These eigenvalues will be between minus one and plus 


is such that this is true, if b is between zero and one. To demonstrate 


this, recall that 





h = Fl F(0), 
-A DUCTU. 
i0 - 
T р 
and hace = Pt) 
Then 1-52 (1-b 
ss 
| NAS 
Then h' = (hy, ho) [1-b2,+(1=b)27. “Hence, 
| | 2 
h +h, = 2(1-b)°, 


which is less than four and is greater than zero. 


Therefore, the vector ont 001 vanishes in the linear model for 


all values of b which are between zero and one. Then 


ЕССЕ) = (aptas) f CO TL nre (1) n 


га (Dro eme DILL ағ ПА me (oT 


h. 


Thus the forecast is unbiased if 


"(D ем а ac 


Τ 


and ea) ΟΕ +h TERRI TUM 


For a specific example see Chapter V. 


C. EXPONENTIAL SMOOTHING COMPARED WITH KALMAN FILTER 

The Exponential Smoothing equation is a special case of the Kalman 
Filter in the steady state. Recall that the matrix F was defined to be 
β μμ Β. κ ҮШ is equivalent to the matrix о” used in the Kalman 
Filter equations, then ве! is identical to the matrix Co derived in 
Chapter II. Thus, in the case of fitting functions consisting of poly- 
nomials, exponentials, and sinusoidals, which can be generated by a 
transition matrix, the Kalman Filter estimate becomes identical with 
that of Exponential Smoothing when Co reaches the steady state. However, 
since the Kalman Filter method does not require that the fitting func- 
tions be generated by means of a transition matrix, it may be applied to 
models other than those listed above. 

The primary difference between the Kalman Filter and Exponential 
Smoothing is the models for which they are valid. The Kalman Filter is 
valid both in a dynamic system as well as in the steady state. The 


Exponential Smoothing technique is strictly valid only in the steady 


state. 


IV. FORECASTING MARINE LOSSES 


In order to see how the Kalman Filter and the Exponential Smoothing 
techniques compared, both methods, along with linear regression, were 
used to forecast U. S. Marine Corps losses in the Republic of Vietnam. 
The actual data used included one hundred and nineteen weeks from 
December 26, 1965, through April 6, 1968. The losses included those 
Marines who either were killed or were wounded and evacuated from 
Vietnam. 

Three models were investigated. These were a growing sinusoidal 


with harmonics with a cycle of sixteen, thirty-two, and fifty-two weeks, 


MEE - a,nta,tt(a,ta,t)sin(pt) + (аа +) со5 (рЕ) 


+ a,sin(2pt) + agcos (2pt) j 
a linear model, 


G€ LA = a, + ant , 


and a constant model, 


The objective was to minimize the mean absolute error where the absolute 
error for the cen forecast was the absolute value of C(t) - C(t). The 
mean error and the estimated variance of the errors were included for 
comparison. 

Both the Kalman Filter and Exponential Smoothing require an estimate 
of the coefficients a(0) in order to begin forecasting. Captain Paul 


William O'Brien, U. S. Marine Corps, had studied the problem of 


Exponential Smoothing in forecasting Marine monthly losses. [O'Brien 1968]. 


23 


Using a least squares line through the monthly losses from March 1965 

to December 1966 and the least squares Fourier fit through the deviations 
about the regression line, Captain O'Brien obtained initial values of the 
coefficients for the growing sinusoidal model. The coefficients of the | 
least squares line through the data were used as the initial values in the 
linear model. These values were scaled by one-fourth and used as the ini- 
tial values for forecasting weekly losses. Thus the growing sinusoidal 


mode | was ιν ized with 
a (0) = (-36.45, 15.675, 62.61, 31.3975, .6325, .6325, 


Б Ио 655 
and the linear model with 


0) = (-36.45, 15.675) . 


A. FORECASTING WITH EXPONENTIAL SMOOTHING 
In each model the value of the weighting factor by which the errors 


were discounted was chosen to be the solution to 
l-b? = 25, 


where n was the number of coefficients in the model. The weighting 
vector h used in the Exponential Smoothing equation was computed by 


determining the steady state matrix F [Brown 1963]. Then 


The transition matrix in each model was determined by inspection, 
rememberina trignometric identities for sin(A+B) and cos(A+B). 
The forecasting procedure involved several steps. These were: 


Step 1; Initialize the vector a(0) and compute f(1). 


24 


Step 2; Compute the forecast C(1) by 


(1) = fl(1)a(0. 


^ 


Step 3; Compute the forecast error 


Step 4: Obtain a new estimate a(1) by 


ad) = L'a(0) + he(1). 


~ ^t 


Step 5; Compute the forecast ο by 


(2) - Ғ(1)а(1). 


-. 


Step 6; Repeat Steps 3 through 5 for each successive forecast 
forecast with the appropriate time arguments. 
The first model investigated was the growing sinusoidal with 


harmonics. For this model the transition matrix is 


0 0 0 0 0 0 0 
1 0 0 0 0 0 0 


0 cosp sinp 0 0 0 0 


L = 0 -sinp cos p 0 0 0 0 | 
0 cosp sinp COST sap 0 0 | 
0 Winep eos pr -sin pWesse pr 0 0 | 
0 0 0 0 O cos ap эт 2 | 
0 0 0 0 0 -sin 2p cos 2p 





where p is 2n divided by the cycle length. Also 


1} = (1, 1, sin p, cos p, sin p, cosgp, sin 2ρ.α TER 


The values of h used were, for the sixteen week cycle 


h? = (0.064709, 0.0011416, 0.025939, 0.12583, 0. 00046034. 


^o 


0.0022076, 04020051, 0-0594527), 


for a thirty-two week cycle 


n! - (0.072061, 0.0012067, 0.056675, 0.12816, 0.00096984, 


- 


0.002 ATL 00497 Ba 


and for a fifty-two week cycle 


h! 


~ 


(0.089418, 0.0013569, 0.10333, 0.13127, 0.0016449, 
0.0020314, 0.063796, 0.029308). 


The results obtained for this model for the three cycles are contained in 
Tables I through III. 
The second model was the linear model. The same procedure was followed 


as outlined above. The transition matrix was 


1 0 
L = 
] | 
and ra) Е (1, 1) 9 
SAA е 


The results from this model are contained in Table IV. 
The last model to which Exponential Smoothing was applied was the 


constant model. With this model 
L = 1], 
μμ η 


and П ж 0275. 
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WEEK 


о о мео ол Боом 


GENERAL EXPONENTIAL SMOOTHING RESULTS 
GROWING SINE WITH HARMONICS 
CYCLE LENGTH 16 WEEKS 


LOSSES 


6] 
64 
84 
23 
87 
116 
МӨ 
108 
167 
405 
89 
85 
539 
208 
80 
116 
14] 
69 
125 
158 
164 
82 
203 
199 
133 
198 
151 
140 
175 
359 
280 
1255 
Cis 
150 
257 
999 
258 
290 
382 
547 
11» 
125 
85 
54 
343 


TABLE I 


ESTIMATED LOSSES 


47 

32 

30 

58 

89 
120 
128 
102 

55 

36 
107 
107 
139 
251 
294 
2/0 
245 
229 
227 
244 
282 
319 
318 
Sy 
321 
304 
315 
302 
252 
352 
405 
404 
087 
491 
348 
307 
402 
411 
418 
426 
45] 
353 
317 
341 
408 


ERROR 


14 
32 
54 
-5 
-2 
-4 
-9 


112 
369 
-18 
-22 
200 
-43 
-214 
-154 
-104 
-96 
-102 
-86 
-118 
-237 
-115 
-198 
-188 
- 106 
-164 
-182 
-157 


-125 
849 
-354 
-341 
-91 
252 
-144 
-121 
-36 
121 
-332 
-228 
-232 
-287 
-0b 


LOSSES. 


118 
96 
80 

262 

131 
85 

177 

106 

435 

145 
9] 

126 

442 
86 

12] 

153 

327 

ІРІ 

220 

336 

280 

237 

202 

300 

273 

389 

646 

732 

359 

400 

253 

247 

249 

248 

536 

249 

256 

254 

35 

208 

269 

203 

290 

608 

297 

260 

189 

166 

262 

216 


TABLE I (CONTINUED) 


ESTIMATED LOSSES 


28 


556 
580 
541 
467 
447 
408 
995 
434 
444 
51 
439 
338 
214 
337 
320 
340 
375 
429 
403 
367 
364 
397 
399 
369 
407 
412 
419 
460 
486 
397 
362 
336 
350 
379 
376 
426 
363 
324 
333 
406 
464 
529 
52] 
495 
517 
439 
360 
317 
321 
371 


ERROR 


-438 
-484 
-46] 
-205 
-316 
-323 
-218 
-328 
-9 
-367 
-348 
-212 
168 
-25] 
-199 
-187 
-46 
-238 
-183 
-3] 
-84 
-120 
-157 
-69 
-134 
-23 
227 
212 
-127 


-109 

-89 
-101 
-125 

160 
-177 
-107 

-70 

-18 
-198 
-195 
-322 
-231 

113 
-160 
-179 
-171 
-15] 

-59 
-155 


WEEK KUSSES 
96 224 
97 304 
98 493 
99 234 

100 262 

10] 466 

102 367 

103 357 

104 414 

105 JUS 

106 243 

10/ 438 

108 208 

109 408 

110 948 

111 2135 

112 1056 

ΜΞ 181 

114 1889 

115 453 

116 563 

117 475 

118 676 

119 982 


Mean Absolute Error 
Mean Error 


Error Variance 


TABLE I (CONTINUED) 


170 
-68.97 
4.42 X 10 


4 


EG 
29 


ESTIMATED LOSSES 


385 
366 
343 
364 
JUS 
338 
440 
500 
520 
305 
444 
ey 
ЭТӘ 
277 
333 
309 
985 
656 
990 
667 
499 
449 
442 
517 


ERROR 


-16] 
-62 
150 

-130 
-61 
128 
-73 

-143 

-106 

-112 

-201 
101 

-]11 
127 
615 
206 
47] 
131 
749 

-214 


26 
234 
65 


TABLE II 


GENERAL EXPONENTIAL SMOOTHING RESULTS 
GROWING SINE WITH HARMONICS 
СҮСВЕ LENGTH 32 WEEKS 


WEEK LOSSES ESTIMATED LOSSES ERROR 
] 61 60 1 
64 64 Р 0 
3 84 63 2] 
4 53 67 -14 
5 8/ 65 22 
6 116 71 39 
7] 119 101 18 
8 108 125 -17 
9 167 145 22 
10 405 177 228 
1] 89 266 1 
12 85 24] -156 
18 339 210 129 
14 208 246 -38 
15 80 23] -151 
16 116 176 -60 
17 141 140 1 
18 163 120 13 
19 125 108 Τη 
20 158 106 52 
2] 164 125 39 
22 82 153 - 
23 203 162 4) 
24 139 212 220 
25 133 240 -107 
26 198 264 -66 
27 15] 301 -150 
28 140 316 -176 
29 1125 32] -146 
30 359 33] 28 
3] 280 386 -106 
32 1253 404 849 
33 23 681 -408 
34 150 619 -469 
35 a 521 2261 
36 92 466 93 
37 258 S05 2208 
38 290 457 -167 
39 382 433 JEG 
40 547 446 101 
4] 119 509 -390 
42 [25 444 2319 
43 85 385 -310 
44 54 346 -292 
45 343 302 4) 


30 


TABLE II (CONTINUED) 


WEEK LOSSES ESTIMATED LOSSES ERROR 
46 118 348 -230 
47 96 324 -228 
48 80 297 -217 
49 262 2/0 -δ 
50 131 300 -169 
5] 85 290 -205 
52 177 ZA -94 
53 106 287 - 181 
54 435 287 148 
95 145 385 -240 
56 9] 389 -298 
57 126 38] -255 
58 442 388 54 
59 86 485 -399 
60 121 464 -343 
61 153 454 -301 
62 327 45] -124 
63 191 492 -301 
64 220 480 -260 
65 336 471 -135 
66 280 486 -206 
67 237 473 -236 
68 202 442 -240 
69 300 401 -101 
70 273 390 -117 
7] 389 573 16 
72 646 30] 255 
73 732 476 256 
74 359 564 -205 
2 400 524 -124 
76 253 496 -243 
17 247 429 -182 
78 249 374 -125 
79 248 334 -86 
80 536 310 226 
8] 249 3/8 -129 
82 256 358 -102 
83 254 347 -93 
84 315 340 -25 
85 208 352 -144 
86 269 331 -62 
8/ 203 331 -128 
88 290 S2 -22 
89 608 520 288 
90 397 415 -58 
SÅ 260 420 - 160 
92 189 394 -205 
23 166 290 -187 
94 262 219 -55 
95 216 32] -105 


3 


WEEK LOSSES 
96 224 
97 304 
98 493 
99 234 

100 262 

10] 466 

102 367 

103 357 

104 414 

105 BUS 

106 243 

107 438 

108 208 

109 404 

110 948 

111 35 

112 1056 

ШЕ 787 

114 1339 

115 453 

116 563 

117 475 

113 6/6 

ΠΠ 586 


Mean Absolute Error 
Mean Error 


Error Variance 


TABLE II (CONTINUED) 
ESTIMATED LOSSES 


161 
-/4.73 
3.929 10 


1 


320 
ο. 
366 
45] 
446 
445 
495 
498 
489 
488 
474 
417 
423 
366 
279 
544 
617 
754 
776 
926 
78] 
678 
564 
53] 


ERROR 


-96 
-25 
127 
-217 
-184 


-128 
-141 
-75 
295 
-23] 


-215 
38 
569 
19] 
439 
39 
563 
-473 
-218 
-203 
ша 
9] 


WEEK 
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GENERAL EXPONENTIAL SMOOTHING RESULTS 
GROWING SINE WITH HARMONICS 
CYCLE LENGTH 52 WEEKS 


LOSSES 


61 
64 
84 
53 
87 
116 
ШЕ 
108 
167 
405 
89 
85 
339 
208 
80 
116 
141 
185 
125 
158 
164 
82 
203 
139 
135 
198 
151 
140 
175 
399 
280 
1253 
273 
150 
257 
559 
258 
290 
382 
547 
117 
125 
95 
54 
343 


TABLE III 


ESTIMATED LOSSES 


33 


63 

73 

79 

87 

83 

88 
102 
114 
12 
147 
235 
210 
208 
260 
269 
238 
220 
211 
200 
187 
185 
185 
161 
175 
169 
163 
177 
176 
174 
165 
244 
2/3 
569 
505 
466 
435 
491 
451 
426 
432 
484 
408 
349 
294 
245 


ERROR 


-2 
-9 


-34 


28 
17 
-6 
44 

258 
-146 
-134 

136 

-52 
-189 
ШТ; 

-79 

-78 

-75 

E 

E 

= 162 
42 

-36 

-36 


-26 
- 36 


176 

36 
980 
-296 
- 385 
- 209 
124 
-233 
- 161 
-44 
115 
- 365 
- 283 
-264 
-240 


Q 
C 





TABLE III (CONTINUED) 


WEEK LOSSES ESTIMATED LOSSES ERROR 
46 118 292 -174 
47 96 275 -179 
48 80 259 -179 
49 262 247 15 
50 151 292 -161 
al 85 296 2511 
De IE 287 ЕО 
53 106 307 220] 
54 435 302 133 
55 145 389 5244 
56 91 376 -285 
57 126 346 2220 
58 442 5227 115 
59 86 397 2311 
60 12] Dor 229] 
6] [55 320 -167 
62 927 302 25 
63 191 335 -144 
64 220 324 -104 
65 336 324 12 
66 280 358 278 
67 231 372 2055 
68 202 579 2:2] 
69 300 366 -66 
70 213 389 -116 
7] 389 403 på 
Да 646 449 197 
73 732 559 173 
74 359 6/2 -313 
75 400 657 -257 
76 293 649 -396 
77 247 596 -349 
78 249 545 -296 
79 248 50] -253 
80 536 462 74 
8] 249 509 -260 
82 256 466 -210 
83 254 43] - 177 
84 315 402 -87 
85 208 394 -186 
86 269 359 -90 
87 203 347 -144 
88 290 320 -30 
89 608 323 285 
90 557 416 -59 
9] 260 422 -162 
92 189 398 -209 
93 166 356 -190 
94 262 318 -5] 
95 216 306 -90 


TABLE III (CONTINUED) 


WEEK LOSSES ESTIMATED LOSSES ERROR 
96 224 287 -63 
97 304 275 29 
98 493 289 204 
99 234 355 E 
100 262 338 -76 
101 466 331 135 
102 367 383 -16 
103 357 398 -4] 
104 414 407 7 
105 393 429 -36 
106 243 440 -197 
107 438 406 32 
108 208 432 -224 
109 404 389 15 
110 948 407 541 
11] 735 576 159 
112 1056 653 403 
113 787 801 -14 
114 1339 840 499 
115 453 1016 -563 
116 563 901 -338 
117 475 825 -350 
118 676 730 -54 
119 582 702 - 180 


Mean Absolute Error = 154 


Mean Error = -67.78 


Error Vi lande = 3.94 X 10" 


14 


m 
σι 


LOSSES 


6] 
68 
84 
56 
87 
116 
119 
108 
167 
405 
89 
85 
η 
208 
80 
116 
14] 
53 
125 
158 
164 
82 
203 
139 
198 
198 
151 
140 
175 
559 
280 
1053 
273 
150 
257 
299 
258 
290 
382 
947 
119 
25 
85 
54 
343 
118 


TABLE IV 


LINEAR MODEL 


ESTIMATED LOSSES 


36 


GENERAL EXPONENTIAL SMOOTHING RESULTS 


11 

42 

65 

88 

94 
109 
128 
142 
148 
170 
257 
227 
201 
258 
260 
222 
203 
196 
187 
177 
179 
182 
159 
178 
73 
166 
181 
77 
ША 
177 
236 
257 
568 
504 
418 
384 
449 
406 
382 
391 
448 
360 
295 
233 
177 
223 


ERROR 


90 
22 
19 
-35 
-7 


-9 
-34 


235 
- 168 
-142 

138 

-50 
-180 
-106 

-62 

-63 

-62 

-19 

-15 
-100 

44 

-39 

-40 

32 

-30 

-37 


182 


996 
-295 
-355 
- 161 

175 
- 191 
-116 


156 
-329 
-235 
-210 
- 179 

166 
-105 


WEEK 


47 
48 
49 
50 
5] 


53 
54 


56 
57 
58 


60 
61 
62 
63 
64 
65 
66 
67 


69 
70 
7] 
73 
74 
75 
76 
71 
78 
79 


8] 
82 


84 
85 
86 
8/ 
88 
89 
90 
9] 
92 
93 
94 
95 
96 


FOSSES 


96 
80 
262 
131 
85 
Т 
106 
435 
145 
91 
126 
442 
86 
| 2 
158 
За: 
[EN 
220 
336 
280 
23/ 
202 
300 
273 
389 
646 
732 
299 
400 
κ. 
247 
249 
248 
536 
249 
256 
254 
Sis 
208 
269 
203 
290 
608 
357 
260 
189 
166 
262 
216 
224 


TABLE IV (CONTINUED) 


ESTIMATED LOSSES 


37 


189 
ve 
128 
162 
148 
124 
134 
120 
210 
190 
159 
146 
232 
190 
168 
161 
209 
204 
209 
248 
261 
ДД 
243 
262 
269 
308 
416 
524 
49] 
477 
422 
ЭЛ 
343 
317 
385 
349 
324 
304 
30/ 
ЯШ! 
2/3 
250 
299 
363 
365 
321 
294 
254 
254 
239 


ERROR 


-93 
-77 
134 
-31 
-63 
53 
-28 
315 
-65 
-99 
-33 
29 
-146 
-69 
-15 
166 
-18 
16 
127 


-24 
-55 
57 


120 
338 
316 
-165 
-9] 
-224 
- 175 
-128 
-95 
219 
- 136 
-93 
-/0 
11 
-99 
-8 
-/0 
40 
349 
-6 
-105 
-148 
-128 


- 38 
-15 


WEEK EOSSES 
97 304 
98 493 
99 234 

100 262 

101 466 

102 367 

103 357 

104 414 

105 393 

106 243 

107 438 

108 208 

109 404 

110 948 

111 735 

112 1056 

113 787 

114 339 

(А15 453 

116 563 

117 475 

118 6/6 

IS 582 


Mean Absolute Error 
Mean Error 


Error Variance 


TABLE IV (CONTINUED) 


128 
-2.2] 
3.607010 


4 


38 


ESTIMATED LOSSES 


231 
249 
324 
296 
286 
340 
S5] 
357 
37/7 
38/ 
348 
377 
330 
353 
536 
609 
761 
794 
984 
857 
793 
716 
8 


ERROR 


73 
244 
-87 
-34 
180 

27 


57 


-144 
90 
-169 


595 
199 
447 


545 
-531 
-294 
-318 

-40 
-136 


TABLE V 


GENERAL EXPONENTIAL SMOOTHING RESULTS 
CONSTANT MODEL 


WEEK LOSSES ESTIMATED LOSSES ERROR 
] 6] 0 6] 
2 64 35 29 
3 84 47 37 
4 53 61 -8 
5 87 58 29 
6 116 67 49 
/ 119 8] 38 
8 108 92 16 
9 16/ 96 7] 

10 405 ІШЕ 290 

11 89 191 -102 

12 85 164 -79 

13 527 144 195 

14 208 194 14 

15 80 197 Hel 

16 116 168 -52 

17 141 155 - № 

18 133 15] -18 

19 125 147 22 

20 158 141 17 

21 164 145 19 

22 82 150 -68 

23 203 133 70 

24 139 151 -12 

25 SS 148 -15 

26 198 144 54 

27 151 158 -7 

28 140 156 - 16 

29 175 152 23 

30 359 158 20] 

280 208 72 

32 1253 226 1027 

33 2/3 483 -210 

34 150 430 -280 

35 257. 360 -103 

36 559 334 225 

ny 258 39] -133 

38 290 397 -6/ 

39 382 341 4] 

40 547 221 196 

41 119 400 -281 

42 125 330 -205 

43 85 279 -194 

44 54 230 -176 

45 343 186 157 
46 118 225 107 


TABLE V (CONTINUED) 


WEEK LOSSES ESTIMATED LOSSES ERROR 
47 96 199 2103 
48 80 173 -93 
49 262 150 112 
50 131 178 -47 
51 85 166 -81 
52 177 146 3] 
53 106 154 -48 
54 435 142 293 
55 145 215 -70 
56 9] 198 -107 
57 126 171 -45 
58 442 160 282 
59 86 230 ЕАД 
60 121 194 -73 
61 153 176 Ко) 
62 ο 170 157 
63 191 209 2189 
64 220 205 15 
65 336 209 127 
66 280 240 40 
67 237 250 -13 
68 202 247 -45 
69 300 236 64 
70 273 252 2] 
7] 389 257 132 
72 646 290 356 
73 732 379 353 
74 359 467 -108 
75 400 440 -40 
76 253 430 2027 
idi 247 386 -139 
78 249 35] -102 
79 248 326 - 78 
80 536 306 230 
8] 249 364 -115 
82 256 335 -79 
83 254 315 EG 
84 315 300 15 
85 208 304 -96 
86 269 280 E 
87 203 277 E24 
88 290 259 3] 
89 608 266 342 
90 357 352 5 
9] 260 353 -93 
2/4 189 330 Е И 
93 166 295 5129 
94 262 262 0 
95 216 262 -46 
96 224 251 297 
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WEEK CDS SES 
oy 304 
98 493 
99 234 

100 262 

101 466 

102 367 

103 397 

104 414 

105 393 

106 243 

107 438 

108 208 

109 404 

110 948 

111 785 

1012 1056 

115 787 

114 1339 

1415 453 

116 563 

117 475 

118 676 

ше 582 


Mean Absolute Error 
Mean Error 


Error Variance 


TABLE V (CONTINUED) 


122 
20.26 
3.48 X 10 


4 


4] 


ESTIMATED LOSSES 


244 
239 
318 
297 
288 
332 
34] 
345 
362 
370 
338 
363 
324 
344 
495 
299 
680 
707 
865 
762 
712 
653 
659 


In this model the initial estimate of a(0) was taken to be zero. Also, 
Since the initial forecasts are biased by E the method of finite 
exponential smoothing [Bessler and Zehna 1968] was used. With this 
procedure, the weight given the forecast error in determining the new 


value of a(t) is h(1-ht)7! 


rather than the constant h. This method 
increased the initial forecasts and as h* vanished, the forecasts 
became identical with those using the constant h throughout. These 


results are in Table V. 


B. FORECASTING WITH THE KALMAN FILTER 
The Kalman Filter requires, in addition to the initial value of 
a(0), the value of c in order to forecast. In the models investigated, 


~ 


the matrix С, was not known. However, since the basic model was 
С = Ва + у 


and in forecasting one time period in advance, D is a scalar, it was 
decided that С, could be estimated as the variance of the forecast error. 
Accordingly, the current estimate of the error variance was used as the 
value of С, in the Kalman Filter equations. 

The initial value of Co was obtained from the equation 

ο, 5 Bier. 

This required at least two forecasts in order to estimate Cy: So that 
a better estimate could be made, Co was initialized after ten time 
periods. The forecasts for these ten periods were obtained by using 
the initialized value of a(Q). 

The forecasting procedure was: 


Step 1; Initialize a(10) = a(0). 


Step 2; Obtain the forecasts for the first ten periods from 
| 

C(t) = b (t)a(10). 

Step 3; Compute the forecast errors, e(t), for these ten 


periods and the mean error, e. 


Step 4; Obtain an estimate of C, from 


m 2 
С = — I [e(k)- e]. 
у ] E 
k=] 
Step 5; Initialize Co by 
Е Ta-l „4-1 
C, (10) - [B C, E 


Step 6; Obtain new forecast CD by 


ο - b (11)a(10). 


Step 7; Compute forecast error e(11) and the new error mean, e. 


Step 8; Obtain new estimate of С, from 


1] 2 
> BREE 
k=] 


= 
ο i 


Step 9; Compute weighting vector K by 


eye παω] 


к = C, (10)b(T1)[C А b 


ED 
= V 2 


Step 10; Obtain new estimate a(11) from 
a dy = Зао еи 
Step 11; Compute new value of C) from 


ΠΣ kb! 


2 (11)c, (10). 
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Step 12; Repeat Steps 6 through 11 with the appropriate 
time arguments for each successive forecast. 

In general the value of С, will be different from "E An exception 
is in the constant model where they are identical. In this model the 
computation of K in Step 9 above involved the old estimate of the error 
variance as the value of Co and С, was the most recent estimate of the 


forecast error variance. 


C. FORECASTING WITH LINEAR REGRESSION 

The method of Least Squares was applied to the data using the 
linear model. The procedure employed was to fit the least squares line 
through the first twelve weekly loss figures. The coefficients of this 
line were used to obtain a forecast of losses for the thirteenth period. 
A new line was fitted to the thirteen observations and the new coeffi- 
cients used for the next forecast. This procedure was continued through 


each succeeding time period. 
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TABLE VI 


KALMAN FILTER RESULTS 
GROWING SINE WITH HARMONICS 
CYCLE LENGTH 16 WEEKS 


WEEK LOSSES ESTIMATED LOSSES ERROR 
] 6] 47 14 
2 64 28 36 
3 84 19 65 
4 53 37 16 
5 87 76 11 
6 116 115 ] 
7. 119 129 -10 
8 108 107 1 
9 167 58 109 

10 405 10 395 

11 89 -7 96 

12 85 384 -299 

18 339 -191 530 

14 208 1146 -938 

15 80 -269 349 

16 116 -590 706 

17 141 189 - 48 

18 183 539 - 406 

19 125 468 - 343 

20 158 285 -127 

2] 164 189 -25 

22 82 100 -% 

28 203 -61 264 

24 139 80 59 

25 133 91 42 

26 198 131 67 

27 151 231 -80 

28 140 255 -115 

29 175 250 -75 

30 359 255 124 

3] 280 278 2 

32 1253 242 101] 

98 273 603 - 330 

34 150 44] -29] 

35 257 272 -15 

36 559 22] 338 

87 258 322 -64 

38 290 275 15 

39 382 236 146 

40 547 226 32] 

4] 119 273 -154 

42 125 19] -66 

43 85 183 -98 

44 54 223 -169 


45 343 288 55 


TABLE VI (CONTINUED) 


WEEK LO SSES ESTIMATED LOSSES ERROR 
46 118 ` 435 ΤΙ 
47 96 453 -357 
48 80 411 2391 
49 262 339 - 77 
50 Б 340 -209 
5] 85 330 -245 
52 177 342 SIGE 
53 106 389 2209 
54 435 387 48 
55 145 422 427] 
56 9] 92] -230 
57 126 197 = 
58 442 118 324 
9] 86 178 ЖО? 
60 12] 180 -59 
61 155 220 -67 
62 327 261 66 
63 191 306 - 115 
64 220 7/2 -55 
65 336 234 102 
66 280 237 43 
67 2] 256 О 
68 202 293 zog 
69 300 332 Е? 
70 СТО 374 = 
71 389 366 2n 
72 646 339 307 
13 732 33] 40] 
74 359 329 30 
76 400 268 132 
76 233 258 -5 
77 247 253 -6 
78 249 265 -16 
79 248 272 224 
80 536 260 276 
8] 249 293 -44 
82 256 265 =o 
83 254 269 -15 
84 315 312 3 
85 208 393 -185 
86 269 448 -179 
87 203 484 -28] 
88 290 464 200: 
89 608 419 189 
90 257 404 547 
91 260 346 Bor 
92 189 30] -112 
93 166 283 207 
94 262 287 -25 
95 216 309 -93 
96 224 301 E 


46 


WEEK LOSSES 
97 304 
98 493 
99 234 

100 262 
101 466 

102 36/ 

103 397 

104 414 

105 393 

106 243 

107 438 

108 208 

109 404 

110 948 

11] 139 

182 1056 

113 787 

114 1359 

115 453 

116 303 

117 475 

118 676 
IM O9 582 


Mean Absolute Error 
Mean Error 


Error Variance 


li 


TABLE VI (CONTINUED) 


166 
11885 
6.28 X 10 


4 


ESTIMATED LOSSES 


279 
257 
290 
308 
953 
459 
510 
Silks 
476 
406 
318 
295 
274 
310 
409 
444 
485 
463 
903 
451 
448 
459 
492 


ERROR 


3] 
236 
-56 
-46 
103 
-92 

-153 
-99 
-83 

-163 
120 
-87 
130 
638 
326 
612 
302 
8/6 
-50 
116 


217 
90 


TABLE VII 


KALMAN FILTER RESULTS 
GROWING SINE WITH HARMONICS 
CYCLE LENGTH 32 WEEKS 


LOSSES ESTIMATED LOSSES ERRCR 


WEEK 
] 61 60 ] 
2 64 63 ] 
3 84 63 2] 
4 93 61 -8 
5 87 62 25 
6 116 69 47 
7 Π5 82 ЭЙ 
8 108 102 6 
9 167 [27 40 
10 405 156 249 
1] | 89 184 -95 
12 85 -75 160 
| 339 79 160 
14 208 155 -947 
15 80 617 -537 
16 116 969 -853 
[ 141 -362 205 
18 52 162 -29 
ІЗ [25 246 - 121 
20 158 183 -25 
2] 164 [19 E 
22 82 117 -35 
23 203 -102 305 
24 1 SK, MNS -76 
25 193 153 -20 
26 198 165 33 
ey 51 345 -]94 
28 140 266 -126 
29 75 227 -52 
30 η gal 68 
31 280 608 -328 
Ва 93 505 748 
E 273 562 -1509 
34 150 1007 -857 
35 297 258 -] 
36 D9 -12 SU 
oy 298 198 60 
38 290 34 258 
38 382 2] 36] 
40 947 144 403 
4] 19 384 - 265 
42 |/25 156 -31 
43 95 50 35 
44 54 -2] T9 
45 343 -67 410 


TABLE VII (CONTINUED) 


LOSSES 


118 
96 
80 

262 

131 
85 

177 

106 

435 

145 
О 

126 

442 
86 

12] 

153 

327 

191 

220 

336 

280 

23 / 

202 

300 

213 

389 

646 

732 

999 

400 

253 

247 

249 

248 

536 

249 

256 

254 

315 

208 

269 

203 

290 

608 

397 

260 

189 

166 

262 

216 


49 


ESTIMATED LOSSES 


102 
5 
47 
24 
89 
85 
ИЗ 

105 

124 

240 

278 

308 

25] 

466 

480 

492 

495 

518 

482 

436 

412 

36/ 

315 

263 

244 

22] 

245 

324 

411 

396 

384 

332 

279 

229 

184 

206 

168 

142 

127 

134 

133 

156 

177 

222 

225 

371 

391 

29 

38] 

3/9 


ERROR 


16 
2] 
33 
238 
42 


104 


311 
-95 
-187 
-182 


- 380 
-359 
-339 
- 168 
-327 
-262 
- 100 
-132 
-130 
-113 

97 

29 
162 
40] 
408 
-52 


-13] 
-85 
-30 

19 
352 
43 
88 
112 
188 
74 
136 
47 
Las 
386 
32 

-11] 

-202 

-225 

-119 

- 163 


WEEK LOSSES 
96 224 
97 304 
98 493 
29 234 

100 262 

101 466 

102 367 

103 357 

104 414 

105 393 

106 243 

107 438 

108 208 

109 404 

110 948 

111 85 

112 1056 

116 787 

114 1339 

ME 453 

116 563 

117 475 

118 676 

119 582 


Mean Absolute Error 
Mean Error 


Error Variance 


18985 
7.64 
9.00 10 


4 


с 


0 


TABLE VII (CONTINUED) 


ESTIMATED LOSSES 


36 1 
340 
312 
358 
334 
322 
347 
306 
365 
382 
393 
374 
384 
351 
348 
419 
144 
500 
504 
965 
504 
460 
414 
404 


ERROR 


-137 
-36 
16] 

-124 
- 72 
144 

20 


49 
11 
- 150 


-176 
53 
600 
3]6 
612 
287 
835 
-112 
59 


262 
178 


TABLE VIII 


KALMAN FILTER RESULTS 
GROWING SINE WITH HARMONICS 
CYCLE LENGTH 52 WEEKS 


WEEK LOSSES ESTIMATED LOSSES ERROR 
1 6] 63 Eo 
2 64 74 -10 
3 84 82 2 
4 53 98 -4] 
5 87 ~45 132 
6 116 -79 195 
7 119 14 105 
8 108 200 97 
9 167 272 -105 

10 405 357 48 

11 89 738 -@о 

12 85 550 466 

NS 339 90 240 

14 208 112 96 

15 80 -78 158 

16 116 -322 438 

17 141 -375 516 

18 155 - 39] 524 

19 125 -569 694 

20 158 -1997 255 

2] 164 8194 -8030 

22 82 7445 -7363 

23 203 5585 -5382 

24 139 3774 -3635 

25 133 2474 -2341 

26 198 1713 -1515 

27 151 1335 - 1184 

28 140 1058 -918 

29 175 850 -675 

30 359 706 -347 

31 280 710 -430 

32 1253 597 656 

33 273 1214 -94] 
34 150 799 -649 

33 257 398 -141 

36 559 192 367 
37 258 242 16 
38 290 65 225 
39 382 -47 429 

40 547 -90 637 
4] 119 -75 194 
42 125 -295 420 
43 85 -537 622 
44 94 -989 1043 


45 343 -3151 3494 


TABLE VIII (CONTINUED) 


WEEK LOSSES ESTIMATED LOSSES ERROR 
46 118 4022 -3904 
47 96 1478 -1382 
48 80 960 -880 
49 262 755 -493 
50 181 803 -672 
5] 85 751 -666 
52 ши 698 -521 
53 106 127 -621 
54 435 708 -273 
95 145 877 - 780 
56 91 836 -745 
57 126 777 -651 
58 442 747 -305 
59 86 865 -779 
60 12] 788 -667 
6] 155 740 -587 
62 327 710 -383 
63 ПЕД 747 -556 
64 220 720 -500 
65 336 703 - 367 
66 280 720 -440 
67 237 710 0713 
68 202 684 -482 
69 300 649 - 349 
70 275 641 -368 
τ] 389 623 - 224 
да 646 631 15 
73 132 689 43 
74 359 750 -39] 
"b 400 724 = 324 
76 253 707 - 454 
77 247 663 -416 
78 249 624 =395 
79 248 590 -342 
80 536 561 -25 
8] 249 583 -334 
82 256 559 - 503 
83 254 540 -286 
84 315 525 2210 
85 208 522 -314 
86 209 504 2235 
87 203 498 - 295 
88 290 481 -191 
89 608 479 129 
90 357 527 2170 
9] 260 527 -267 
92 189 509 -320 
93 166 478 2212 
94 262 443 - 181 


95 216 423 -207 


TABLE VIII (CONTINUED) 


WEEK LOSSES ESTIMATED LOSSES ERROR 
96 224 495 -171 
97 304 368 -64 
98 493 351 142 
515) 234 859 -125 

100 262 328 -66 

101 466 302 164 

102 367 303 64 

103 357 291 66 

104 414 280 134 

105 3» 278 115 

106 243 276 -33 

107 438 257 181 

108 208 2/0 -62 

109 404 429 149 

110 948 271 677 

111 75 358 377 

112 1056 414 642 

113 787 906 28] 

114 1339 560 779 

IRIS 453 675 -222 

116 963 677 -114 

117 475 690 -215 

118 676 689 -13 

119 582 710 -128 


Mean Absolute Error = 65] 


Mean Error - 359.16 


Error Variance E07. 07 10? 


TABLE IX 


KALMAN FILTER RESULTS 
LINEAR MODEL 


WEEK LOSSES ESTIMATED LOSSES ERROR 
1 61 22] 
| 3 22 
4 29 26 27 
5 87 4] 46 
6 116 57 59 
7 119 73 46 
8 108 88 20 
9 167 104 63 
10 405 120 285 
1] 89 135 -46 
12 85 122 397 
13 339 127 212 
14 208 185 23 
15 80 209 - 129 
16 116 203 287 
17 141 203 -62 
18 133 206 E 
19 125 206 28] 
20 158 204 -46 
21 164 208 -44 
22 82 210 -128 
28 203 200 З 
24 139 209 -70 
25 5 207 -74 
26 198 203 E 
27 15] 210 -59 
28 140 209 -69 
29 175 206 = 
30 359 208 15] 
31 280 234 46 
22 1253 247 1006 
33 273 286 TB 
34 150 294 -144 
35 297 298 E 
36 559 305 254 
37 258 324 -66 
38 290 330 -40 
39 382 337 45 
40 547 348 199 
4] 119 366 -247 
42 125 364 2090 
43 85 36] 2975 
44 54 357 2505 


45 343 35] _8 


WEEK 


46 
47 
48 
49 
50 
9] 
52 
53 
54 
39 
56 
57 


59 
60 
6] 
62 
63 
64 
65 
66 
67 
68 
69 
70 
7] 
72 


74 
75 
76 
27 
78 
ща 
80 
8] 
82 
83 
84 
85 
86 
87 
88 
89 
90 
2 
92 
93 
94 
98 


LOSSES 


118 
96 
80 

262 

131 
85 

177 

106 

435 

145 
91 

126 

442 
86 

121 

155 

327 

19] 

220 

336 

280 

237 

202 

300 

23 

389 

646 

732 

359 

400 

253 

247 

249 

248 

536 

249 

256 

254 

315 

208 

269 

203 

290 

608 

55 7 

260 

189 

166 

262 

f 


TABLE IX (CONTINUED) 


ESTIMATED LOSSES 


358 
355 
350 
344 
347 
343 
537 
355 
329 
340 
336 
329 
324 
355 
328 
322 
318 
SES 
321 
320 
325 
327 
326 
324 
3 7 
328 
334 
393 
373 
378 
383 
382 
380 
378 
377 
388 
386 
384 
382 
383 
380 
3/0 
375 
305 
388 
390 
388 
384 
3/9 
3 


ERROR 


-240 
-259 
-2/0 

-82 
-216 
-258 
- 160 
-229 

106 
- 19 
-245 
-203 

118 
-249 
- 207 
-169 


-132 
- 101 
16 
-45 
-90 
-124 
-24 
-54 
61 
ἽΝ 
379 
-15 
22 
-130 
-135 
-131 
- 130 
189 
-139 
-130 
-130 
-6/ 
- 175 
- 111 
-176 
-85 
235 
-3] 
- 130 
- 199 
-218 
-117 


WEEK LOSSE? 
96 224 
97 304 
98 493 
29 234 

100 262 

101 466 

102 367 

103 387 

104 414 

105 398 

106 243 

107 438 

108 208 

109 404 

110 948 

un 735 

112 1056 

113 787 

114 1359 

MIS 453 

116 563 

117 475 

118 6/6 

ІШ? 582 


Mean Absolute Error 
Mean Error 


Error Variance 


a 
— 


TABLE IX (CONTINUED) 


141 
- 15:88 
4.18 АЩО 


4 


56 


ESTIMATED LOSSES 


374 
3/2 
Sie 
380 
За 
376 
382 
385 
387 
391 
394 
391 
396 
392 
395 
416 
429 
450 
462 
485 
488 
494 
497 
505 


ERROR 


-150 
-6δ 
12] 

-146 

-115 

90 
-15 
-28 

27 


-151 


-188 
Ше 
558 
3119 
627 
ser 
ӘЛІ 
-32 
75 
-19 
179 
77 


TABLE X 


KALMAN FILTER RESULTS 
CONSTANT MODEL 


WEEK LOSSES ESTIMATED LOSSES ERROR 
] 61 0 6] 
2 64 0 64 
3 84 0 84 
4 53 ] Bip 
5 87 29 58 
6 116 97 39 
7 119 80 39 
8 108 103 5 
0 167 103 64 

10 405 163 242 

1] 89 164 6 

12 85 117 - 

1 339 94 245 

14 208 136 72 

15 80 208 -128 

16 116 208 -92 

|. 14] 149 -8 

18 133 142 a 

19 125 137 E 

20 158 13] 27 

2] 164 157 7 

22 82 158 = 

23 203 153 50 

28 139 201 -62 

25 183 198 265 

26 198 165 30 
27 151 198 -47 
28 140 196 -56 
29 175 171 4 
30 359 175 184 
3] 280 176 104 
32 1255 260 993 
33 273 267 6 
34 150 273 -123 
35 257 264 -7 
36 559 258 30] 
37 258 271 23 
38 290 259 3] 
39 382 287 05 
40 547 305 242 
4] 119 317 -198 
42 [25 247 510? 
43 85 164 -79 
44 54 pies -58 
45 343 7 266 


97 


TABLE X (CONTINUED) 


WEEK LOSSES ESTIMATED LOSSES ERROR 
46 118 12 -5 
47 96 118 202 
48 80 110 280 
49 62 96 166 
50 131 141 -10 
JE 85 133 -48 
52 ІШ? 120 57 
93 106 176 =70 
54 435 174 261 
55 145 228 285 
56 9] 166 -75 
57 126 125 ] 
58 442 126 Зб 
59 86 132 - 10 
60 121 90 а 
6] 23 120 33 
62 327 140 187 
63 ІР! 141 50 
64 220 191 29 
65 336 194 142 
66 280 198 92 
67 207 270 333 
68 202 264 -62 
69 300 242 58 
70 273 299 -26 
7] 389 298 91 
72 646 360 286 
13 732 37] 361 
74 359 5065 -146 
15 400 401 =1 
76 238 400 A 
77 247 39] -144 
78 249 317 -68 
79 248 266 EE 
80 536 252 284 
8] 249 268 -19 
82 256 250 6 
83 254 255 E 
84 315 204 6] 
85 208 305 -97 
86 269 304 -35 
87 203 277 -74 
88 290 254 36 
89 608 290 318 
90 257 290 67 
9] 260 356 -96 
92 189 398 -164 
95 166 302 236 
94 262 223 39 
95 216 262 -46 


98 


WEEK 


96 

97 

98 

YE 
100 
101 
102 
188 
104 
105 
106 
107 
108 
109 
110 
111 
ІШЕ 
ШӘ 
114 
15 
116 
117 
118 
Це 


224 
304 
493 
234 
262 
466 
36/ 
357 
414 
393 
243 
438 
208 
404 
948 
735 
1056 
787 
1839 
453 
563 
475 
676 
582 


Mean Absolute Error 


Mean Error 


Error Variance 


LOSSES 


ESTIMATED LOSSES 


109 
41.33 
ОЗЕ ТО 


5 
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TABLE X (CONTINUED 


262 
24] 
296 
302 
254 
261 
267 
355 
SED 
400 
398 
ЗЯ 
438 
437 
406 
415 
660 
813 
788 
798 
581 
563 
548 
637 


ERROR 
- 38 


197 
-68 


205 
100 


58 
-/ 
-155 
47 
-230 
-33 
542 
320 
396 
-26 
55] 
-345 
-18 
-88 
128 
-55 


TABLE XI 


LINEAR REGRESSION RESULTS 


WEEK LOSSES ESTIMATED LOSSES ERROR 
13 239 196 143 
14 208 252 -44 
15 80 256 -176 
16 116 224 -108 
17 14] 208 ~67 
18 133 201 -68 
18 125 ІШЕ -68 
20 158 194 225 
2] 168 184 220 
2: 82 185 108 
23 203 170 33 
24 139 179 -40 
25 133 176 -43 
26 198 [72 26 
27 151 178 227 
28 140 177 297 
29 175 174 ] 
30 359 177 182 
3] 280 203 77 
За 11253 217 1036 
33 213 350 277 
34 150 351 -201 
35 257 337 -80 
36 559 336 223 
37 258 369 -111 
38 290 366 276 
39 382 366 16 
40 547 376 17] 
4] 119 402 -283 
42 125 383 -258 
43 85 366 -28] 
44 54 347 -293 
45 343 327 16 
46 118 334 -216 
47 96 320 -224 
48 80 306 -226 
49 262 29] -29 
50 13] 92 -161 
51 95 283 -198 
52 177 271 -94 
53 106 266 -160 
54 435 257 178 
55 145 272 -127 
56 91 265 -174 
57 6 255 -129 


60 


WEEK 


58 
39 
60 
61 


63 
64 
65 
66 
6/ 
68 
69 
70 
7] 
ща 
13 
74 
15 
76 
II 
78 
79 
80 
81 
82 
83 
84 
85 
86 
87 
88 
89 
90 
9] 
D 
93 
94 
95 
96 
97 
98 
99 
100 
101 
102 
105 


LOSSES 


442 

86 
12] 
153 
327 
191 
220 
336 
280 


а 


202 
300 
273 
389 
646 
732 
359 
400 
293 
247 
249 
248 
536 
249 
256 
254 
ΙΙ; 
208 
269 
203 
290 
608 
357 
260 
189 
166 
262 
216 
224 
304 
493 
234 
262 
466 
367 
357 


TABLE XI (CONTINUED) 


ESTIMATED LOSSES 


6] 


248 
263 
253 
246 
24] 

248 
246 
246 
253 
256 
256 
255 
259 
26] 

2/0 
293 
319 
323 
930 
329 
Sen 
ΕΕ 
324 
337 
336 
334 
338 
334 
33 
51510 
за 
327 
342 
345 
343 
997 
334 
333 
330 
327 
328 
337 
334 
333 
340 
343 


ERROR 


194 
-177 
-132 

-93 


-5/ 
-26 
90 
27 
-19 
-54 
45 
14 
128 
376 
439 
40 
77 
-77 
-82 
-/8 
-/8 
22 
-88 
-80 
-80 
-18 
-126 
-62 
oc] 
-3/ 
281 
[5 
-85 
-154 
-173 
-/2 
-117 
-106 
-23 
165 
- 103 
-/2 
133 
27 
14 


TABLE XI (CONTINUED) 


WEEK BOSSES ESTIMATED LOSSES ERROR 

104 414 346 68 

105 295 350 43 

106 243 354 -111 

10/ 438 352 86 

108 208 397 -149 

109 404 323 3] | 
110 948 SE 59! | 
ТШ > 380 398 
112 1056 395 661 

ІШЕ 787 421 366 

114 1339 437 902 

115 453 472 -19 

116 563 474 89 

117 475 480 -5 

118 6/6 484 192 

119 582 493 89 


Mean Absolute Error = 139 


m 


Mean Error 


Error Variance = 4.53 X 10 
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Figure 1 
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EXPONENTIAL SMOOTHING (LINEAR FORECASTS) 


Figure 2 
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Figure 3 
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24 48 72 96 120 
KALMAN FILTER (CONSTANT) FORECASTS WEEKS 
Figure 4 
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24 48 72 96 120 
EXPONENTIAL SMOOTHING (CONSTANT) FORECAST WEEKS 


Figure 5 


6/ 


X 10 


12 


10 


FORECAST 
LOSSES 


24 48 72 
LINEAR REGRESSION FORECASTS 


Figure 6 
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96 


WEEKS 


120 





V. CONCLUSIONS 


In the case of the linear model, it was shown in Chapter III that 


E[C(t)] = (ала в (п-в e (OY! in 


ms 


аз cren AT 1 a Js 


for large t and when the actual data was such that 


E[C(t)] = a, + apt. 


Then,  E[C(t)] = a, + aot + ao. 


This showed that when the assumed model was correct, the forecast was 
biased by the slope of the data line. In the case of the problem of 
forecasting the Marine losses, the least squares line through the one 
hundred and nineteen points had slope 3.41. This suggested that Exponen- 
tial should have produced a forecast which is 3.41 more than the actual 
data, or an error of -3.41. The actual mean error was -2.21. Also, the 
Exponential Smoothing forecast would be above the linear regression fore- 
cast when the losses were increasing and below when they were decreasing. 
The results bear this out, although not around the time periods when the 
actual losses changed directions. For example, the losses increased 
during the time periods thirty-seven through forty. The Exponential 
Smoothing forecast with the linear model exceeded the linear regression 
forecasts during these periods. The losses then decreased during the 
periods forty-one through forty-four. The Exponential Smoothing forecast 
was above the linear regression forecast for the forty-first week but it 


then was below the linear regression forecasts through forty-four. 
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STANDARD DEVIATION OF ERRORS VS. MEAN ERROR 


Figure 7 


The results obtained with the Kalman Filter and the growing sine model 
were generally comparable with the results with Exponential Smoothing. 
This is evidenced in Figure 7. The exception was the sine model with a 
cycle length of fifty-two weeks. The Kalman Filter forecast for this 
model fluctuated badly for the initial forecasts. However, the last 
seventy weeks of the data produced a mean absolute error of -321.71 and 

a mean error of -208.71. The last forty weeks produced a mean absolute 
error of 219 and a mean error -24.08. When these results are compared 
with the overall mean absolute error of 651 and the mean error of -359.16, 
it appears that the forecasts were improving as the effects of the initial 
forecasts were dampening out. Thus, the model whould not be totally dis- 
missed, especially since the Exponential Smoothing results with the cycle 
length of fifty-two weeks did not vary greatly from the results when a 
cycle length of sixteen and of thirty-two weeks was used. The Kalman 
Filter technique assumes that the forecast error variance is known. For 
the problem investigated, that variance was estimated, which was a variation 
on the Kalman Filter. If the true variance had been known, the results of 
all the models with which the Kalman Filter was used might have been im- 
proved. 

With the models used in Gans SN Marine losses, the simpler the 
model, the better the results. Captain O'Brien found that the linear model 
produced better results than the more complicated growing sine model for 
the prediction of monthly Marine losses [O'Brien 1968]. The same was true 
for the weekly forecasts. In addition, the even er constant model 
produced even better results. However, if the actual cycle period had been 
used rather than an estimate, the results of the Exponential Smoothing with 
the growing sine model might have been greatly improved. The period length 


can be determined as described in R.G. Brown's book[Brown 1963, p. 66-75]. 


7) 


In general, the estimated mean of the forecast errors was not zero. 
For example, the Kalman Filter with the constant model, on the average, 
forecasted forty-one fewer losses than actually occurred. With this 


information, a better policy would have been 


where u is the mean forecast error, С! (+) is the estimated losses and 

C(t) is the forecasted losses from the method used. The variance of the 
forecast error will not be changed by this modification since C(t) would 
still be computed using C(t-1) - C(t-1). However, the mean error, where 


the error is 


will be zero. 

Finally, the results indicated, as evidenced in Figure 7, that even 
when an estimate of the forecast error covariance matrix С, is used, the 
Kalman Filter results do not differ greatly from those of Exponential 
Smoothing. Since, except for the constant model, the Kalman Filter is 
more complicated than either the Exponential Smoothing or the Linear 
Regression methods, either of these latter methods should be used when 
С, is unknown or when no computer facilities are available. When Cy is 
known and a computer is available, the Kalman Filter technique might 
provide better results than either of the other two methods. 

Under the criterion chosen, the minimum mean absolute error, the 
Kalman Filter method using the constant model was the best forecasting 
method. This method could be easily used for forecasting Marine losses 
in Vietnam. The calculations for this constant model are relatively 


simple and can be accomplished with a calculator without the aid of a 


computer. The calculation of the estimate of the forecast error variance 
is the most complicated portion of the constant model. The mean absolute 
error and the variance of the masm. achieved with this model 
compare favorably with the data mean of 294.6 and variance of the data of 
b. 2 = 10°. Using this model would be more beneficial than merely using 
the data mean as the forecast of the losses for the next time period. 


Because of the difference between the two variances, one would be more 


certain of a minimum error with the Kalman Filter constant model. 
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GENERAL EXPONENTIAL SMOOTHING PROGRAM 


This program computes successive weekly forecasts for the model 


E[c(t)] = a, + а Е. 
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